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1. Introduction 

In this note, we study the moments of the average of the Ramanujan sum. To be precise, 
we want to evaluate the sums 

k 



n<y v q<x 

where k is a positive integer, x and y are large real numbers, and 
(1.1) c g (n) = Y e(-an/q) = ^ d ^(l/ d ) 

l<a<q d\(q,n) 
(a,q)=l 

is the Ramanujan sum. Our interest in this question stems from an old attempt [TJ to apply 
Fourier techniques to Diophantine approximations of reals by sums of rational numbers. At 
the time, we were surprised that information on the asymptotic behavior of Ck(x, y) appeared 
to be missing from the literature. Of course, it is easy to evaluate Ck(x, y) in some cases. For 
example, elementary arguments based on the second formula in (11.11) yield the asymptotic 
formulas 

C 1 (x,y) = y + 0(x 2 ^) 

and 

C 2 (x,y) = ^ + 0(x^), 

for any fixed e > 0. However, these bounds (as well as other simple things we have tried) are 
only of interest when y > x 2+e . In our application to Diophantine approximations, on the 
other hand, we were interested in C2(x,y) when x 1+e < y < x 2+e . Some numerical experi- 
mentation suggested that in that range the order of C 2 (x, y) is still yx 2 , though the coefficient 
2^2) seemed somewhat off when y is close to x 2 . In the present note, we prove, among other 
things, that those empirical observations are true. Our first result is an asymptotic formula 
for Ci(x,y). 

Theorem 1.1. Let x be a large real number and y > x. Then 

Ci{x,y) = y- -— + O (xy 1/3 logx + x 3 y^). 

The proof of Theorem 11.11 uses exponential sum estimates and can be easily modified 
to yield different estimates for the remainder. The remainder claimed in the statement of 
the theorem is the result of applying one of the simplest exponential sum bounds — the 
exponential pair (|, |) from van der Corput's method [2]. More sophisticated exponential 
pairs will result in replacing the error term xy 1 ^ 3 by different ones, which will be sharper for 



certain choices of the relative sizes of x and y. The interested reader will have no trouble 
obtaining such improvements by appropriate modifications of (I2.5P below. Since our main 
focus is on maximizing the "support" of the result, we prefer the simple version above. 

Theorem 11.11 demonstrates that when x = o(y), the Ramanujan sum c q (n) is o(l) on 
average over q < x and n < y. Our next theorem shows that this is no longer the case when 
the average C%(x,y) is replaced by the mean square sum C 2 (x,y). 

Theorem 1.2. Let x be a large real number, y > x, and B > 10 be fixed. 

(i) If y > x 2 (logx) B , then 

C 2 (x, y) = ^j + 0{yx\\ogx) w {x-V 2 + xV 1 )). 

(ii) If x(\ogx) 2B <y< a; 2 (log a;)- 8 , then 

C 2 (x,y) = J^(l+2*(«)) + 0(yx 2 {\ogx) w (x- 1 ' 2 + (y/x)- 1 / 2 )), 

where u = \og(yx~ 2 ), and k{u) is defined by (13 .201) below and satisfies the inequalities 

k(u)> -0.4, k(u) < exp ( - |w| 3/5 " e ), 

for any fixed e > 0. In particular, k(u) = o(l) as \u\ — > oo. 

The hypothesis B > 10 and the powers of the logarithms in the ranges for y in the 
statement of the theorem are chosen to ensure that the error terms are of a smaller order 
of magnitude than the main term. Also, the lower bound for k(u) is far from sharp and has 
been chosen so to simplify the arguments in §H In fact, by increasing the amount of numeric 
calculations in §Hby an order of magnitude or two, it can be shown that n(u) > —0.3. We 
chose not to pursue such a sharper bound here for the sake of clarity. However, for the benefit 
of the reader, we should mention that the more elaborate numerical computations that we 
carried suggest that k(u) is decreasing when u < u$ and increasing when u > u$, where 
Uq ~ 1.63; its minimum value is approximately k(1.63) ~ —0.2943. It also seems that k(u) 
is always negative. 

2. The first moment 

Using (11.11) . we have 



n<y q<x d\q n<y dk<x 

d\n d\n 

Interchanging order of summation, we obtain 

(2.i) c x {x, y ) = dmJ2 1 = E ^(*) 



dk<x n<y dk<x 

d\n 



dk<x dk<x dk<x 

Ci,!(x, y) - C lj2 (x, y) - C lj3 (x, y), say. 
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Here, for a real t, [t] is the integral part of t and ip{t) = t — [t] — | is the saw-tooth function. 
It is easy to see that 



m<x k\m 



(2.2) 
Further, 

(2.3) c h2 (x, y) = \Y, Mk) = \Y1 E d 

dk<x k<x d<x/k 



k<x v 7 



X 



4 k 2 

k<x 



4C(2) 



+ 0(x \ogx) 



It remains to estimate C\^(x, y). By partial summation, 

oo 

(2.4) C 1 , 3 (x,y)«^^iV,sup ^(y/d) 



k<x j=0 



del 



where Nj = Nj^ = (x/k)2~ : > and the supremum is over all subintervals I of (Nj,2Nj]. We 
remark that the sum over j is, in fact, finite and has 0(log x) terms. We now use the following 
lemma, which is a special case of [21 Lemma 4.3] (see [2] or (3] for the definition of exponential 
pairs) . 

Lemma 2.1. Suppose that (x, A) is an exponential pair and that 1 is a subinterval of(N, 2N] . 
Then 

We apply Lemma 12.11 with the exponential pair (x, A) = (|,|) to the sum over d on the 
right side of (12.41) . This yields 

oo 

(2-5) C 1>3 (x, y) « )T ('V + NjV~ l ) 

k<x j=0 

« ((*/% 1/3 + V 1 ) 

k<x 

< (xy 1/3 logx + x 3 y~ 1 ). 
Theorem 11.11 now follows easily from (I2.1I) - (I2.5I) . 



3. The second moment 



Recall the identity 



(3.1) 



E 

9=1 



c q {n) = ai- s (n) 
q s C(s) 



(Ms) > 1), 



where £(s) denotes the Riemann zeta-function and a z {n) = J2d\ n d z - We note that without 
loss of generality, we may assume that x,y G Z + ~. Then, using ( 13. ip and the truncated 
Perron formula [3J Corollary 5.3], we get 

V- / x 1 f a+iT <r\-s{n) x s , „ , , 

Here, T is a real parameter at our disposal, a > 1 + (logx) -1 is a real number, and the 
remainder Ei(x,n) satisfies 

9=1 ^ x/2<g<2x ^ 1 

If we assume further that a < 1 + 3(logx) _1 , we can show that 

(3.3) Ei(x,n) <t: (x/T)a (n) log x. 

Applying ( 13. 2 p and (13. 3p with a = acj = 1 + j(logx) -1 , j = 1,2, we obtain 

2 i /-qi+iT i>a.2+iT 



^ c g (n) ) = , ., 2 / / F(si,s 2 ,n)ds 2 dsi + E 2 (x,n), 

„ ,. / Joi-iT Ja 2 -iT 

o"i_ sl (n)o"i_ S2 (n) x sl+S2 



9<x 

where 

i?Sl ' S2,Tl/ C(si)C(s 2 ) SiS 2 

and 

E 2 {x,n) <C (x 2 /T)a (n) 2 logx. 
Summing (13 ,4p over n, we deduce 

(3.5) C 2 [x,y) = —— i \ / ds 2 dsi + — — 

(2m) 2 J ai _ iT J a2 - iT C(si)C(s2j sis 2 V T 

where L = \og(xy) and 

G(si,s 2 ;y) = ^0-i_ Sl (n)0i_ S2 (n). 

n<y 

We now recall Ramanujan's identity 

>^ o- a (ra)(T b (ra) _ C(gXtf ~ a )(( s ~ b )C( s -a-b) 
1 j ^ n s C(2s-a-6) 

n=l v ' 

valid for Re(s) > max{l, 1+Re(a), 1+Re(6), l+Re(a+6)}. From (13. 6p . by another application 
of Perron's formula, we get 

I ra+3iT w fyL A \ 

(3.7) G(si,s 2 ;y) = — H(s u s 2 ,w)—dw + 0[ — ), 

2tt? J a _ 3iT w \ T J 

where a = 1 + 3(logy) _1 , y G Z + ~, and 

C(w)C(w + gi - l)C(w + g 2 - l)C(w + si + g 2 - 2) 

H{si,s 2 ,w) = — ■ ■ — . 

£(2w + si + s 2 - 2) 

Let r(a, 0, T) denote the contour consisting of the line segments [a—iT, (3—iT], [ft—iT, 0+iT] 
and [(3 + iT, a + iT]. We now move the integration in (13. 7p to T(a, ~, 3T); let us denote the 



respective integrals Ji, I2 and I3. When Re(sj) = ctj, the integrals l\ and J3 over the two 
horizontal line segments are bounded above by 

(3.8) J 1 ,/ 3 «C(«i + « 2 -l)L 4 / T l - 2 °y°da<^L\y l l 2 + yT- 1 ). 

J 1/2 

Here, we have used the standard convexity bound (see [U eqn. (5.1.4)] for a slightly weaker 
version) 

(3.9) \((a + ft) I « (|*| + 2)^)/ 2 logflf | +2) (1 < a < 1). 
Furthermore, by Holder's inequality, 

(3.10) h^y^Ci^ + ^-l) J] M 4 (|+j(«i-l) + A;(a 2 -l),5T), 

0<j',fc<l 



where for | < er < 1, 



M 4 4 (a,T)= [ T \((<j + it)\ 

J-T 



1 + 1*1 

Appealing to the fourth-moment estimates for £(s) jH §7.5 & §7.6], we deduce from (13.101) 
that 

(3.11) / 2 «y 1 / 2 L 6 . 

Combining (13.71) . (13. 8p and ( 13. lip , we obtain 

4 

(3.12) G( Sl , s 2] y) = Y, R M> *2l y) + 0{L\y 1 ' 2 + yT' 1 )), 

where Rj(si, s 2 ; y) are the four residues of the integrand in (13. 7p : 

3/C(*i)C(s2)C(si + s 2 - 1) 



-Ri(si,s 2 ;y) 
■Ba(si, s 2 ;y) 
i2 3 (si,s 2 ;y) 
^4(^1, s 2 ;y) 



^1 + ^) 

j / 2 -^C(2-si)C(i-^i + ^)C(^) 

(2- Sl )C(2-si + s 2 ) 

j / 2 -X(2-s 2 )C(i + ^i-^)C(^i) 

(2-s 2 )C(2 + Sl -s 2 ) 
y 3-i-*, C (3 _ Sl _ S2 ) C ( 2 _ S2 ) C ( 2 - Sl 



(3 - si - s 2 )C(4 - sx - s 2 ) 
Substituting (13.12)) into the right side of (I3.5p . we get 

4 

(3.13) C 2 (x,y) = ^C 2j (a;,y) + O( 2/ x 2 L 10 (y- 1 / 2 + T- 1 )), 

3=1 

where 

1 r i-HT n^ Rjjsuwy) ^** J 

C 2 , j {x,y) = 7 —— 2 - / / ds 2 dsi- 

(2ttz) j y Ql _ iT J a2 _ iT C(Sl)C(S2) «iS 2 

We now proceed to evaluate the integrals C 2 j(x, y), starting with C 2 ,i(x, y). We first move 
the integration over s 2 to the contour T(a 2 , \ ,T)\ let us denote the integrals over the three 
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line segments J\ t \, Ji j2 and J\$. The contributions from the integrals over the two horizontal 
line segments are bounded using (13. 9ft : 



< xyL 2 (xT- 1 + 1). 



Furthermore, using estimates for the second moment of ((s) in the critical strip, we find that 
the contribution from the integral over the line Re(s2) = \ is 



^1,2 < / / 777 i -in/1, i T\ dtidt 2 



3/2 / / ISV^l 2 

-ri-T (l*i| + l)(l*2|-- I 

-2T /-T 



dtdu 

T (\t\ + l)(\t - u\ + 1) 



«^ 3 / 2 y ic( ai -j+<«)i y 

/2T j 
IC(«i-| + m)|^-« 2/ x 3 / 2 L 3 . 
-2T M + 1 



Hence, accounting for the residue at s 2 = 2 — si, we have 

^ 2 + 0(yx 2 L 3 (a;- 1 / 2 + T- 1 )). 



2C(2) 

To estimate 6*2,4(0;, ?/), we move the integration over S2 to the contour T(a 2 , (3,T), where 
j3 = | — «i; let us denote the respective integrals J^i, and 74,3. Similarly to the estimation 
of Jij, we have 

Ja,u As « yx 2 L 5 (T^ + (x/yT) 1 / 2 ), J 4 , 2 « y l/2 x 5/2 L\ 
Since the integrand is holomorphic in the region between the two contours, we deduce that 

(3.15) C 2A (x,y) «yx 2 L 5 ((x/ 2 /) 1 / 2 + T- 1 ). 
Similarly, by moving the integration to r(a 2 ; ^T), we find that 

(3.16) C 2 , 3 (x,y) ^yx'L'Ux/y^ + T- 1 ), 
and 

(3.17) C 2;2 (x,y) = ^K(x,y) + 0(yx 2 L 5 ((x/y)^ + T- 1 )), 
where 

1 p+^ C(2 _ s )( yx -2)i- 



K(x,y) = / — — -ds. 

The first term on the right side of (I3.17P arises from the pole of R 2 (si, s 2 ;y) at s 2 = s±. 
Thus, by fl3~T3j) - (l3~T7j) . 

(3.18) C 2 (x,y) = |^(| + K(x,y)) + O (yx 2 L 10 (x^ 2 + (x/yf 2 + T" 1 )). 
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It remains to evaluate the integral K(x,y). To that end, we argue differently depending on 
the size of yx~ 2 . 

Case 1: yx~ 2 > L B for some fixed B > 10. Then we move the integration to the contour 
r(ai,/3,T), where (3 = 2 — L _1 ; let us denote the respective integrals Ki(x,y), K 2 (x,y) and 
K 3 (x, y). Then 

ir 1 (x, 2/ ),K 3 (x,i/)«r- 2 , 

and 

K 2 {x, y) « y- l x 2 L J |C(2 - - it)\ ™ « fVi. 
Therefore, in this case, on choosing T = x, we obtain 

(3.19) C 2 (x,y) = ^ ) +0{yx 2 L w (x^ 2 +x 2 y- 1 )). 

This completes the proof of part (i) of the theorem. 

Case 2: yx~ 2 < L B . Then we move the integration to the contour r(«i,l,T). Let us 
denote the respective integrals K^(x, y), K 5 (x, y) and K$(x, y). Then, by the bounds for £(s) 
in |H §6.19], 

K 4 (x,y),K 6 (:r,y)<T- 2 , 

and 

K 5 (x,y) = n(u) + 0(T- 2 ), 
where u = log(yx~ 2 ) and k(u) is the Fourier integral 

(3.20) .(u)^l£f(U)e-S, /(.) = «i^i| Ir _± r _ y . 
We choose T = x and get the asymptotic formula 

(3.21) C 2 (x,y) = ^-(\ + «(«)) + O^L^x/v) 1 '* + x^ 2 )). 

Thus, to complete the proof of part (ii) of the theorem, we need only establish the desired 
properties of k(u). We examine k(u) in the next section. 

4. Estimation of k(u) 

We now study the function k(u). Suppose first that \u\ — > oo. For a fixed e > and a real 
t, we define 

5(t)=a(log(|t|+2))- 2 / 3 -. 
By the standard bounds for £(s) near the edge of the critical strip [U p. 135], 

1/(^ + ^)1 « (^+^1)3 when |a| < 5(t). 
Applying Cauchy's integral formula for the circle jt, \s — it\ = 8(t), we deduce that 



I /«(«)! 



fc! /" /(*) 



2vri y (s - it) k+1 



ds 
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k\8(t) 



-k-2 



Hence, a /c-fold integration by parts gives 



K[U) < 



\U\ 



f ik \it)e~ ltu dt 



< 



k\C: 



-fc-2 



(logt) (2/3+ e )(fc+2) r 3 dt _ 



Using the inequality 



we deduce that 



(logt) a t _d (it < / :"r 

2 JO 



~ ..-^r^ + i) 



2 a+l 



k\rrik\ 



£,(5/3+e)fc 



2 m *(C e M) fe (C e |M|) fc ' 
where m fc = [(2/3 + e)(A; + 2)]. Finally, we set k = [C^ | ^ | 1/(5/3+e) J , and we get 

- l«|3/5-^ 



< exp ( - |m| 3/5 e ). 



Next, we focus on the case when \u\ is bounded. The integral k(u) is differentiable, and 



WM\ < — 



\t\dt 1 

2tt (i + t 2 ) 3 / 2 - 2- 



Hence, 
(4.1) 



-3 



it — u 



Since the integrand is bounded above by \t\ , we have 

r-60 

-itu 

J V^je 

-60 

where |i? u | < |. Let /Cq(m) denote the integral on the right side of (14. 2p . By partial integration, 



(4.2) 



K(U) 



i r 

^ J -60 



dt + $ u 10~ 



|«o(w)| < 



27T U 



|/(60i)| + |/(-60i)|+ / |/'(it)|dt < 



60 



60 



4.2 

2tt\u\ 



< 0.67|w 



-i 



Combining the latter inequality and ( 14. 2p . we deduce that |k(w)| when |w| > 1.7. It remains 
to estimate k(u) when \u\ < 1.7. In view of ( 14. ip and ( 14. 2p . it suffices to tabulate «o(w) for 
a set of suitably chosen values ui,U2, ■ ■ ■ , with |-Uj| < 1.7. To this end, a quick calculation 
using Mathematica yields the following table. 





K ( Uj ) 


covered range 


-1 


-0.0277 


-1.7 < m < -0.3 


0.3 


-0.0984 


-0.3 < u < 0.9 


1.1 


-0.2586 


0.9 < it < 1.3 


1.5 


-0.2925 


1.3 < u < 1.7 



Each row of this table lists a value of Ko(uj) and a range of the variable u, where we can 
infer the bound k(u) > —0.4 from ( 14.11) . ( 14. 2 p and the listed value of K (uj). For example, 
when 0.9 < u < 1.3, we have 

k{u) > «(1.1) - 0.5|w - 1.1| > Ko(l.l) - 0.1001 > -0.36. 



Clearly, this completes the proof of part (ii) of Theorem 11.21 
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